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Abstract Multiple solutions exist in many experimental situations when several inter- 
fering amplitudes are summed to fit experimentally measured distributions, such as cross 
sections, mass spectra, and/or the angular distributions. We show a few examples where 
multiple solutions are found, but only one solution is reported in the publications. Since 
there is no standard rule for choosing one among the solutions as the physics one, we 
propose a simple rule that agrees with what has been adopted in previous literatures: the 
solution corresponding to the minimal magnitudes of the amplitudes must be the physical 
solution. We suggest test this rule in the future analyses. 
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1 Introduction 


In quantum mechanics, a physics observable is proportional to the squared modulus 
of the amplitude. In case of more than one amplitude contributing to a process, they 
are summed to obtain the total amplitude, and thus one generally has contributions from 
interference terms to the physics observable. It is simple to predict a physics observable 
when the amplitudes and the relative phases between them are known, in which case there 
are no ambiguities. 

However, in many circumstances, the experimental quantities are measured, and from 
these one extracts information on the amplitudes. As there is a square operation on the 
summed amplitudes in calculating the observable, one expects that mathematically, there 
can be multiple solutions. In these cases one faces a common problem: which solution is 


the physics one. 
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In this paper, we present a few examples where two-solutions are reported in fitting 
with the coherent sum of two amplitudes; and revisit a few cases where only one solution 
is reported and where we find other distinct solutions from fits with two amplitudes; we 
also discuss more complicated situations where more than two amplitudes are summed 
to fit the data. Finally, we propose a conjecture on how to choose the physics solution 


among multiple solutions. 


2 Examples with two-solutions reported 


A few recent examples where multiple solutions are reported are studies of the so- 
called Y states that are produced via initial state radiation by the Belle experiment. 
Figure 1 shows the invariant mass distributions of ntn J/1» and a*a~w(2S') measured 
in Belle |1, 2], together with a fit with two coherent resonant terms and an incoherent 
background term. Table 1 shows the fit results, including the Y (4008) and Y (4260) 
Breit-Wigner (BW) terms from the n*m J/w mode, and the Y (4360) and Y (4660) BW 
terms from the »*« v(2S) mode. It should be noted that in both channels there are 
two solutions with exactly the same goodness-of-the-fit (x?/ndf = 81/78 for ntm J/v 
mode [1] and y?/ndf = 4.7/3 for x * v v (2S) mode [2]), with exactly the same masses and 
widths for the resonances but with very different couplings to e*e^ pair (Tete-). Instead 
of choosing one of the two solutions, Belle reported both in their publications; however, 
the Particle Data Group (PDG) frequently uses only one of the solutions in averaging 
with results from other experiments [3, 4]. Such a treatment of data is rather capricious, 
since the solutions picked up randomly from different experiments may have distinctive 
features and, thus, should not be averaged together. 

Another example is the study of the decay dynamics of n — ynn mode by BES [5]. 
When the c *7- invariant mass distribution is fitted with coherent sum of the p resonance 
and a contact term, it is found that there are two solutions of equal goodness-of-the-fit. 
One solution corresponds to constructive interference between the two amplitudes while 
the other to destructive interference. However, in previous analyses only one solution [6] 
was reported, the one that corresponds to the solution with constructive interference in 
the BES analysis [5]. 


3 Examples with one-solution reported 


Since the examples described in the previous section are basically fits to an s-dependent 
distribution with two coherent amplitudes, one would suspect that there would be, in 
general, two solutions in all such circumstances. There are many of these types of fits 


in the literature, especially in low energy ete~ annihilation experiments, where the cross 
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Figure 1: The rtr J/w (upper) and ata v(25) (lower) invariant mass distributions and 
the best fit with two coherent resonances together with a background term. The data are 
from Belle [1, 2]. 


section of e*e^ — hadrons is parameterized as the coherent sum of the amplitudes of the 
vector mesons. Below, we show two typical cases where only one solution is reported, and 


we have redone the fit to obtain the other solution from the data. 


3.1 Branching fraction of ọ — wr” 


The most precise data on e*e^ — wn near the ¢ resonance were reported by KLOE [7] 
for both the w > «x 3? and w — yr? decay modes. In the KLOE paper [7], the cross 


section as a function of the center-of-mass energy, ys, is parameterized as 


a(s) = Onr(v/'8) ` à <i ee 


where On; (1/8) = co +.0'(4/s — Mg) is the bare cross section for the non-resonant process, 
parameterized as a linear function of \/s; Z is the interference parameter, while My, Tg, 
and D, = M, » — s —iM3T, are the mass, the width, and the inverse propagator of the @ 
meson, respectively. 

We take the KLOE data for w — -*- 5? from Table I of Ref. [7] and fit with the 
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Table 1: Fit results of the ata J/v and «v *« wv(2S) invariant mass spectra |1, 2]. The 
first errors are statistical and the second systematic. M, Tot, and B-TI';«,- are the mass 
(in MeV), total width (in MeV), product of the branching fraction to hadronic mode and 
the e*e- partial width (in eV), respectively. $ is the relative phase between the two 


resonances (in degrees). 


Parameters Solution I Solution II 
M(Y (4008)) 4008 + 40*114 
Tot (Y (4008)) 226 + 44 + 87 
B-T.+--(Y(4008)) 5.041476) 12.4 = 2.47148 
M (Y (4260)) 4247 + 1213 
T tot (Y (4260)) 108+ 19+ 10 
B.T.«-(Y(4260) 6.01.2447 20.6 2.3191 
Q 12 + 29*7. —111z 7*2 
M(Y (4360)) 4361+99 
Tot (Y (4360)) 74+1510 
B -Tee (Y(4360)) 1044£1.7+15 11.8 1.8 1.4 
M(Y (4660)) 4664+115 
Piot (Y (4660)) 4815-3 
B-T.+--(Y(4660)) 3.0+0.9+0.3 7.61.8 0.8 
$ 39+30+22  —79 +1720 


parametrization given above. In our fit, the Born-order cross section is calculated as 
Ovis/Orad, Where ovi, is the visible cross section and ó;44 is the radiative correction factor; 
only statistical errors are considered in the y? construction. Table 2 shows the results 
from our fit, where two solutions are found with the same fit quality (x?/ndf = 4.3/13). 
We can see that the parameters from Solution I are very close to those listed in Table II 
of Ref. [7], the slight difference is due to the simplified procedure for the Born-order 
cross section calculation in our fit. It can be seen that Solution II has a much larger 
resonant amplitude, and the resulting branching fraction of the isospin violating process 
$ — wT? is at per mille level, about two orders of magnitude higher than Solution I, 
namely, the solution reported in the original work [7]. Figure 2 shows the fit results 
and the contribution of each component in the fit. Similarly, we made a fit to the cross 
sections measured with w + y7°, there are also two solutions obtained from the fit and 


these results are in good agreement with those from w — ata s. 
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Table 2: Results from fits to the ete~ — wr? cross sections measured with w — ata~7°. 


Parameter Solution I Solution II 

co [nb] 7.88 + 0.04 7.88 + 0.08 
R(Z) 0.106 +0.004 0.106 + 0.006 
SCA) —0.103 + 0.003  —1.90 + 0.006 


o'[nb/MeV]  0.0642-0.002 — 0.064 + 0.006 
B(ó— wm?) 461x105 7.62 x 1073 
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Figure 2: Fit to the ete” > wr 
vs. Data points are shown with dots with error bars. Left plot is for solution I: at 
ys = 1.019 GeV, from top to bottom, the solid curves are continuum, best fit, resonance, 
and interference terms. Right plot is for Solution II: at /s = 1.019 GeV, from top to 
bottom, the solid curves are resonance, continuum, and best fit, the interference term is 


not shown. 


3.2 F, and the p-w mixing 


Even before the discovery of the 7 lepton, Tsai calculated the branching fractions of 
such a heavy lepton decaying into vector final states such as T^ — «c «7v, by applying 


the conserved vector current (CVC) hypothesis to e*e^ — ntr data measured in e*e^ 
CVC 


annihilation experiments [8]. This calculation, referred to as By’S,,-,0,,, has been tested 
ever since the discovery of the 7 lepton. As the precision of both the ete~ annihilation 
and the 7 decay experiments improved significantly in the past three decades, the test has 


reached the precision level of better than 1%, and a discrepancy at the 2c level between 
CVC 


the measured 7^ — v 1?v, branching fraction B,-.,,-,o,, and BEYS 0, is observed [9]. 
Many theoretical efforts have been made to explain this discrepancy, including a better 


understanding of the isospin violation correction. 
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One source of isospin breaking is the effect of p-w mixing in the e*e^ annihilation data 
eve 


tT 07 Wp, 


ete” — ata data in the vector meson dominance (VMD) model (including p and its 


that is absent in 7 decays. The mixing effect on B is estimated by fitting the 
excited states and w) and subtracting the w contribution by setting its amplitude to zero. 
This has been carried out in many previous analyses [9, 10, 11, 12], but in all cases only 
one solution is reported. As a heuristic example, we only use the CMD2 data and follow 
the fit described in the CMD2 publication [13], namely, with the GS-parametrization [14] 


of the pion form factor F}: 


BW$5(s) [1 + 63: P.(s)] + 8BWEP(s) 


Fos S 
pu t 
where 
2 . 
BWE(s) = Utd Tvem) 
mg — s+ f(s) — imyT'v(s) 
2 
m 
P, = = . 
(s) m2 — s — iml, 


The definitions of all the quantities can be found in Refs. [11, 12, 14]. Here, the complex 
numbers ô and 6, as well as the mass and width of the p, are fit parameters. The masses 
and widths of w and p’ are fixed to their PDG values [3, 4], and the phase of £ is fixed at 
180°. 

We find two solutions with equally good quality (x?/ndf = 21/25) of fit, as shown in 


Fig. 3 and Table 3. It is clear that Solution I is in good agreement with the results reported 
CVC 


by CMD2 analysis [13], and the resulting p-w mixing corrections to B% 5,9, and a, 


(muon magnetic anomaly) are in agreement with the results in Refs. [9, 10]. However, 


the relative w and p strengths are different in the other solution, as are the corrections to 
peve 


T 07 Wp, 


As the amplitude of the w term in Solution II is more than an order of magnitude 


and a,, as shown in Table 3. 


larger than that in Solution I, we would expect a branching fraction of w — ntr much 
larger than that reported in previous analyses: about 7096 in Solution II versus 1.596 in 
Solution I. It should be emphasized that B(w — ntr?) depends on the measurement 
of w — atm, w — yr", and so on, since they all are derived from ete~ annihilation 
experiment, and there are no independent measurements of B(w — nrn?) [3, 4]. The 
well known B(w — ntra’) z 90% is a result of taking Solution I of the fit to F;. 
Solution II seems to be in contradiction with expectations that the isospin breaking decay 


w — ntn should be much smaller than the isospin conserving decay w > mtron’. 
However, it is also the case that we have come to accept a small  — ntr branching 


fraction because the existence of the second solution in the fit has never been considered. 


IF, 
IF, 
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Figure 3: Fit to the e*e~ — mtx form factors below s = 1 GeV? measured at CMD2 [11, 
12]. Left: along the data points (dots with error bars) is the best fit. The / — 1 part of 
the fit is shown for the two solutions. Right: comparison of the p-w interference part in 


two solutions. 


4 Examples with more than two amplitudes 


The examples shown above are all described by the sum of two amplitudes with an 
unknown relative phase in the fit. In the case of more than two amplitudes contributing to 
the distribution, the number of solutions is 2”~', where n is the number of free amplitudes. 

In the fit to the ete” — mtr cross sections [15], incoherent and coherent sums 
of two amplitudes (the $(1680) and Y (2175)) are considered, while in estimates of the 
significance of a third resonance ( X (2400)), Belle also tested the possibility of the coherent 
sum of three amplitudes. This procedure is repeated with combined BaBar [16, 17] and 
Belle data [15] in Ref. [18]. In this latter case, four solutions are found with the same 
masses and widthes for the resonances, but different coupling constants. 

The same phenomenon was also shown in the fit to combined BaBar [19] and Belle [2] 
data on ete~ — «*« v(2S) [20], where four solutions were found when fitting with the 
coherent sum of Y (4360), Y (4660), and Y (4260) amplitudes. 


5 Partial wave analysis 


All the examples discussed so far are one-dimensional problems. In the case of fitting 
a multi-dimensional distribution, the same multi-solution problem also exists. The partial 
wave analysis (PWA) and the Dalitz plot analysis techniques used widely in hadron physics 


potentially have such a problem. 


Table 3: Results from fits to the e*e^ — ntn” form factors measured at CMD2 [11, 12]. 


AB™™ng and Aq S are the corrections to BOYS and a,,, respectively, due to p-w 


T —m- m9 Vr 


mixing. 
Parameter Solution I Solution II Davier [9] 
mp [MeV] 715.9 + 0.5 - 
T, [MeV] 146.0 + 0.8 - 
ló| [x1073] 1.62+0.06 — 21.97 + 0.04 - 
ós [°] 10.1+1.4 — 86.56 0.17 - 
|6| 0.086 0.004 - 
AgBmixing [6] — —0.03 +0.01 +0.04+0.01 | —0.01 + 0.01 
Agmixne [190719] 42.5402 — 416-0. | +2.80 0.19 


In a recent BES analysis of the (25) — «*«- J/V process [21], the 1 *«^ system is 
fitted with contributions from c, a contact term, and a possible tensor amplitude f2. It is 
curious that the two very similar models shown in Figure 4 of Ref. [21] give very different 
strengths for the c term. The aforementioned expositions led us to suspect that two differ- 
ent solutions are found in each of the two models, one with constructive interference and 
the other with strong destructive interference. Fortunately, according to the study above, 
we found that all different solutions result in an identical line shape for the resonance (the 
same mass and width), and, thus, the resonant parameters are unaffected. However, since 
the fraction of each component can be very different for different solutions, the branching 
fraction for each component obtained from these analyses can not be compared directly 
if they are from different solutions. 


A typical partial wave analysis has many amplitudes and all of them are added coher- 
ently to fit the data distributions. In principle, there are multiple solutions, although we 
cannot judge how many-fold the ambiguity may be. In some circumstances, when more 
than one solution is found, a simple average has been taken to be the best estimate of the 
parameters. Such a treatment is obviously non-physical, since any one of the solutions is 
a good description of the data but the average is generally not. In some other cases, the 
multiple solution problem is just overlooked, which makes the comparison between exper- 
iments not meaningful, since one compares two different solutions from two experiments 
which are different by definition. Therefore, some care must be taken in comparing the 
branching fractions obtained from fitting multi-dimensional distributions with the coher- 
ent sum of several amplitudes. Unless they are from the same solution, the comparison 


is meaningless. 


6 Minimal amplitude conjecture 


In the preceding sections, we examined a few cases where amplitudes are extracted 
from experimentally observed distributions, and in all these cases more than one solution 
was found. We also discussed the similar problem in PWA. All these indicate that this 
problem actually is common to many analyses. However, this problem is often overlooked 
and, more often than not, only one solution is reported with no indication of any reason 
for neglecting the others. Moreover, as we believe the physics must correspond to only one 
of these solutions, there is no existing rule or solid physics argument on how one should 
select one solution from all the possibilities. 

However, from the existing examples and our observation, we notice that the solution 
with the smallest modulus of the amplitude always agrees with our expectation, for ex- 
ample the isospin violating processes ¢ — wr? and w — mtr should be small if w is an 
isospin zero dominant state. So it may not be very surprising that the physics solution 
of any system corresponds to the one where the modula of the amplitudes take the min- 
imal values among all the solutions. We call this selection rule “the minimal amplitude 
conjecture". 

Such a rule can be comprehended readily from a simple “economic” principle. One 
would expect that a solution in which the physics observable is produced via two or more 
very large amplitudes with strong destructive interference is less economic than the one 
via small amplitudes but with constructive interference. In the latter case, we always find 
a maximum constructive interference among amplitudes. With this minimal amplitude 
conjecture, we can determine the unique physics solution of all the aforementioned ex- 
amples. This will make theoretical analyses applicable to cases where many solutions are 


extracted from experimental data. 


7 Conclusion 


In this paper, we put forth two important problems: multiple solutions in fitting 
experimental data and the selection of the physics solution among many possible solutions. 
As a matter of fact, more and more cases with multiple solutions appear in recent data 
analyses and even more may have been overlooked previously. Here, we strongly suggest 
that all possible solutions be found out and reported in the future analyses, in order to 
avoid a babel of arguments and misleading theoretical deductions. 

In theoretical analysis, when confronting with multiple solutions, one has to make 
a choice. We propose a “minimal amplitude conjecture”: the solution with minimal 
modula of the amplitudes to be the physics one. Such a rule ensures a unique solution 


be singled out from multiple solutions. This conjecture could further be tested with more 


experimental data analyses. If more and more facts are found to support it, it could be 


promoted to a principle. 


Note added: We were reminded by the Journal referee of a preprint arXiv:0710.5627 [22], 
where the author studied mathematically the ambiguity in determining the resonant pa- 
rameters in fitting to the cross sections with the coherent sum of two or more Breit-Wigner 
functions, and obtained analytical relations between different solutions for some special 
forms of the amplitudes. 
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